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ABSTRACT: In this paper, we study the Peristaltic 
transport of magnetohydrodynamic (MHD) Jeffrey fluid in 
a non-uniform porous channel with the influence of slip, 
wall properties and heat transfer under the assumptions 
of long wavelength and low Reynolds number. The 
analytical expressions for the stream function, velocity 
and temperature are obtained. The results for velocity, 
stream function and temperature obtained in the analysis 
are discussed through graphs. It is noticed that the 
velocity and temperature decrease with increasing Jeffrey 
parameter \ . Further it is observed that the size of the 

trapped bolus decreases with increasing i, . 
Keywords: 
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INTRODUCTION 

Peristalsis is an important mechanism for mixing 
and transporting fluids, which is generated by a 
progressive wave of contraction or expansion 
moving on the wall of the tube. This mechanism is 
found in the swallowing of food through esophagus, 
chyme motion in the gastro-intestinal tracts, 
movement of ovum in the fallopian tube and many 
other glandular ducts in a living body. The 
mechanism of peristaltic transport has been 
exploited for industrial applications like sanitary 
fluid transport, roller and finger pumps, blood 
pumps in heart lung machine and transport of 
sensitive or corrosive fluids where the contact of 
the fluid with the machinery parts is prohibited. 

The inertia free peristaltic flow with long 
wavelength analysis was given by Shapiro et al., 
[1]. The early development on mathematical 



modeling and experimental fluid mechanics of 
peristaltic flow was given by Jaffrin and Shapiro [2]. 
A theoretical study of peristaltic transport of two- 
layered power-law fluids is made by Usha and Rao 
[3]. Kavitha et al., [4] studied the peristaltic flow of a 
Williamson fluid in an asymmetric channel through 
porous medium. Many researchers have 
contributed to the study of peristaltic transport 
under the effect of magnetic field and porous 
channel [5-10]. 

Peristaltic transport in non-uniform ducts is 
considerable interest as many channels in 
engineering and physiological problems are known 
to be of non-uniform cross-section. Srivastava et 
al., [11] and Srivastava and Srivastava [12] studied 
peristaltic transport of Newtonian and non- 
Newtonian fluids in non-uniform geometries. 
Radhakrishnamacharya and Radhakrishna Murthy 
[13] studied the interaction between peristalsis and 
heat transfer for the motion of a viscous 
incompressible fluid in a two-dimensional non- 
uniform channel. Mekheimer [14] studied the 
peristaltic flow of blood (obeying couple stress 
model) under the effect of magnetic field in non- 
uniform channels. He observed that the pressure 
rise for a couple stress fluid is greater than that for 
a Newtonian fluid. Also the pressure rise for 
uniform geometry is much smaller than that for 
non-uniform geometry. Hariharana [15] investigates 
the peristaltic transport of non-Newtonian fluid, 
modeled as power law and Bingham fluid, in a 
diverging tube with different wall wave forms. 

Mittra and Prasad [16] analyzed the peristaltic motion 
of Newtonian fluid by considering the influence of the 
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viscoelastic behaviour of walls. They assumed that the 
driving mechanism is in the form of a sinusoidal wave of 
moderate amplitude imposed on the flexible walls of the 
channel. Dynamic boundary conditions were proposed 
for the fluid motion due to the symmetric motion of the 
flexible walls which were assumed to be either thin 
elastic plates or membranes. Radhakrishnamacharya 
and Srinivasulu [17] studied the influence of wall 
property on peristaltic transport with heat transfer. Sobh 
[18] introduced slip effects on couple stress fluid. 
Ramana Kumari and Radhakrishnamacharya [19] 
investigated the effect of slip on peristaltic transport in 
an inclined channel with wall effects. The influence of 
slip, wall properties on MHD peristaltic transport of a 
Newtonian fluid with heat transfer and porous medium 
have been investigated by Srinivas and Kothandapani 
[20]. 

Among several non-Newtonian models proposed 
for physiological fluids, Jeffrey model is significant 
because Newtonian fluid model can be deduced 
from this as a special case by taking 1, = 0. Further 

it is speculated that the physiological fluids such as 
blood exhibit Newtonian and non-Newtonian 
behaviors during circulation in a living body. The 
Jeffrey model is relatively simpler linear model 
using time derivatives instead of convective 
derivatives for example the Oldroyd-B model [21]. 
Kothandapani and Srinivas [22] studied the 
peristaltic transport of a Jeffery fluid under the 
effect of magnetic field in an asymmetric channel. 
More recently, Vajravelu et al., [23] studied the 
influence of heat transfer on peristaltic transport of 
a Jeffrey fluid in a vertical porous stratum. 

Motivated by these studies, in the present 
investigation we study the Peristaltic transport of 
magnetohydrodynamic Jeffrey fluid in a non 
uniform porous channel with the influence of slip, 
wall properties and heat transfer under the 
assumptions of long wavelength and low Reynolds 
number. The closed form of solution for the 
velocity, stream function and temperature are 
obtained. The effects of different physical 
parameters on the velocity, stream function and 
temperature obtained in the analysis are discussed 
through graphs. 



BASIC EQUATIONS 

The constitutive equations for an incompressible 
Jeffrey fluid are 

T=-pT+Y (1) 

J = — ( f + A, y) 

i+V ; 

where T and J are Cauchy stress tensor and 
extra stress tensor respectively, p is the pressure, 

/ is the identity tensor, X x is the ratio of relaxation 
to retardation times, X 2 is the retardation time, f is 

shear rate and dots over the quantities indicate 
differentiation with respect to time. 

MATHEMATICAL FORMULATION 

We consider the motion of an incompressible 
electrically conducting Jeffrey fluid in a two- 
dimensional non-uniform porous medium channel 
induced by sinusoidal waves propagating with 
constant speed 'c'. The walls of the channel are 
assumed to be flexible and are taken as a 
stretched membrane. The wall deformation h(x,t) 
due to the infinite train of peristaltic waves is 
represented by 

y = h(x,t) = d (x) + asm — (x-ct) (3) 

X 

where d(x) = d + mx, m«l, a is the 

amplitude, X is the wavelength, c/ is the mean half 
width of the channel, m is the dimensional non- 
uniformity of the channel. 



c 


X 














i<7 h(x.t) 

• • • v '. 
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Fig. 1 Physical Model 
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The equations governing the flow of an 
incompressible Jeffrey fluid in a porous medium 
under the influence of a magnetic field are 

:0 



du dv 
dx dy 



V 



du du du 

YU h V 

dt dx dy 



J 



dx 



as„ ds r 



dx dy 



(4) 



(5) 



dv dv dv 

YU h V 

dt dx dy 



-GB n u u 

0 k 

dy dx 



dS, 



dy 



v(6) 



and the energy equation [24] is 



CP 



dT dT dT 

YU Y V- 

dt dx dy 





( d 2 T | d 2 T^ 


= K 

) 


v dx 1 dy 2 / 


du 






S — + S 
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5 »VM 



du du 

1 

dx dy 

(7) 



where 

S„ = 



2// 



1+4 



1 + 4 
2// 



1+4 



1+4 



1+4 



1+4 



d d^ 

u 1- v — 

dx dy 



d d 

u 1- v — 

dx dy 



d d^ 
u ^-v— 

dx dy j 



du 
dx 

du dv 
— + — 
dy dx 

dv_ 

dy 



Here «,vare the velocity components along x 
and y directions respectively, p is the density, is 
the coefficient of viscosity of the fluid, p is the 
pressure, a is the electrical conductivity of the 
fluid, B 0 is the intensity of the magnetic field acting 
along the y-axis and the induced magnetic field is 
assumed to be negligible, k is the permeability of 
the porous medium, £ is the specific heat at 
constant volume, v is the kinematic viscosity of 



the fluid, K is the thermal conductivity of the fluid 
and T is the temperature of the fluid. 

The governing equations of motion of the flexible 
wall may be expressed as 

L t {h) = p- P{) (8) 

where L* is an operator, which is used to 
represent the motion of stretched membrane with 
viscosity damping forces such that 

d 1 d 2 d 
L= - T d? +mi d? + C dt 
Continuity of stress at 
momentum equation yields 



(9) 

y = ±h and using x 



dx dx dx 



dy 



V 



du du 

YU 

dt dx 



du 



— oBlu u 

k 



_, du 



u = +h at y = ±h = ± 



d + m'x + a sin 



In 
T 



dy 

(10) 
[x—ci) 
(11) 



T = T„ on y: 



-h 



T = T X on y = h (12) 
Herer is the elastic tension in the membrane, 
m is the mass per unit area, C*is the coefficient of 
viscous damping, p Q \s the pressure on the outside 
surface of the wall due to the tension in the 
muscles and /z, is the dimensional slip parameter. 

It is assumed that p 0 =0 

We introduce the stream function ^"such that 

dw dw 
u = — v = — . 

dy dx 

and the following non-dimensional quantities are 

, x , y , w , ct 
x =—, y =—, y/ = — , t =—, 
Pi d cd Pi 

, d T-T 0 , h , d 1 

S = — 5, 6 = -, h = — , p = p 

juc T t -T Q d cXp 



(13) 
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The non-dimensional governing equations after 
dropping the primes, we get 



RS 



d 2 y/ 3 y/ d 2 y/ 3 y/ d 2 y/ 
dtdy dy dxdy dx 3y 2 



dy 



RS 3 



d 2 y/^dy/d 2 y/ dy/ d 2 y/ 
dtdx dy dx 2 dx dxdy 



dp „dS n 
— - + 8 — — 
dx dx 

dy/ \ dy/ 
dy Da dy 
(14) 

dy dx 
S 2 dy/ 



RS 



3<9 dy/ d6 _ dy/dd 

dt dy dx dx dy 



1 

Pr 





' 3 y/ 






^3x3_y 


V 



V 



dy Da dx 
(15) 
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(17) 



dx dy 
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3x 3 
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(18) 



Further, it is assumed that the streamline value is 
zero at y = 0 . i.e. ^(0) = 0 
(19) 

<9 = 0 on y = -h and (9 = 1 on y = h (20) 

where 

a „ <i 

4 1 



are 



geometric parameters, 
^_cdp_ j s ^ ne Reynolds number, M = J— B 0 d is 



the 



Hartmann 



I// 

number, 



-Q/ 3 



are the non-dimensional elasticity parameters, 



Pr 



number, m 



is the Prandtl number, £c = — ; 

k 

— is the Darcy 

d 2 

Xm 



is the Eckert number, Da 



is the non-uniform parameter 



and /? is the Knudsen number (Slip parameter). 

EXACT ANALYTICAL SOLUTION 

Under the assumptions of long wavelength ( 
8 «1) and low Reynolds number, from equations 
(14)-(18), we get 



0: 



0 = 



+ - 



1 



dp 

dx ' (1 + 4) dy 3 

dp 

dy 



.<*V_ M 2 ^ 1 dy/ 

dy Da dy 



Equation (22) implies p * p(y) 



0 = ^? + 



Ec (d 2 y/^ 2 



Vxdy 2 1 + 4 I 3y 



(21) 
(22) 



(23) 



(h) 



Elimination of pressure from equations (21) and 
(22), yields 

(24) 



aV A^ = 0 



3/ 



dy 2 
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where N = J(1+A [ ) 



M 2 + — 
Da 



Equation (18) gives 



dy 3 dy ' 



dx dxdt dxdt 



(h) 



(25) 



The closed form solution for stream function and 
velocity from equation (24) using the boundary 
conditions (17), (19) and (25) are given by 



¥- 



Sex 3 



N 



[E^ + E 2 ) cos 2n{x -t)- sin 2n{x - 1) 



x 



sinhAfy 



N ( cosh Nh + N/3 sinh Mi ) 



-y 



(26) 



Sett 3 



u =- 



N 2 



[E { +E 2 )cos2ft(x-t)-^- sin 2n(x - 1) 



x 



cosh Ny 



(cosh M + AT/? sinh M) 



— 1 



(27) 

Substituting equation (26) into equation (23) 
subject to the boundary condition (20), we get the 
temperature as 

Brl\ (2N 2 y 2 - cosh 2Ny + cosh 2Nh -2N 2 h 2 ) ( y + h) 



0 



8(1 +4)( cosh Nh + N/3 sinh Nh) 2 



2h 
(28) 



where 

Uk 3 



N 



Y~ sin 2n(x -t)-(E 1 + E 2 ) cos 2n{x - 1) 



and Br= Ec Pr is the Brinkman number. 



The coefficient of heat transfer at the wall is given 

byz = h x -e y 



(m + 2ft£cos2x(x-t)) 



8h (cosh Nh + N/3 sinh Nh) 



/i (l + /i 1 ) L ' ^ 4Nly ~ 2N Sinh 2Ny ) + 4 ( cosh Nh + N fi sinh Nh f 



(29) 



RESULTS AND DISCUSSION 

The equation (27) gives the expression for the 
velocity in terms of y. Velocity profiles are plotted in 
figures from (2) to (8) to study the effects of 
different parameters such as the Jeffrey parameter 
\, the Darcy number Da, the slip parameter , 
the amplitude ratio e , the non-uniform parameter 
m , the Hartmann number M, the wall tension E l , 
the mass characterizing parameter E 2 and the 
damping nature of the wall E 3 . Fig. 2 is drawn to 
study the effect of \ on the velocity distribution u. 

We observe that the increase in \ decreases the 
velocity. 



Fig. 3 is plotted to study the effect of Da on the 
velocity. We observe that the velocity increases 
with increasing Da. Further, the permeability 
parameter causes to strengthen the fluid slip at the 
wall. Fig. 4 and Fig. 5 shows that an increase in /? 
and e results in the increase of velocity 
distribution. In Fig. 6 we find that the velocity for a 
divergent channel (m>0) is higher compared to 
its value for a uniform channel (m = 0), whereas it 
is lower for a convergent channel (m<0). From 
Fig. 7 we see that the velocity decreases with 
increase of M. From Fig. 8 we notice that the 
velocity increases with increasing E, and E 2 and it 
decreases with increasing E 3 . 

The equation (28) gives the expression for the 
temperature in terms of y. Temperature profiles are 
plotted in figures from (9) to (16) to study the 
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effects of the physical parameters of the problem. 
Fig. 9 is drawn to study the effect of \ on the 
temperature distribution 6 . It is observed that the 
temperature decreases with increasing \ . Figures 
1 0 and 1 1 are plotted to study the effect of Da and 
M on the temperature. We observe that the 
temperature increases with increasing Da and 
decreases with increasing M. Figures 12 and 13 
shows that the temperature decreases by 
increasing f3 and it increases with increasing e . 

Fig. 14 illustrates the effects of m on the 
temperature distribution. We notice that the 
amplitude of the temperature is large in case of 
divergent channel compared with uniform and 
convergent channels. Fig. 15 is plotted to study the 
effect of Brinkman number Br on the temperature 
distribution. We notice that the temperature 
increases with an increase in Br. Fig. 16 shows 
that the temperature increases with increasing E, 
and E 2 and it decreases with increasing E 3 . 

The equation (29) gives the expression for the 
coefficient of heat transfer at the wall. Figures from 
(17) to (23) are plotted to observe the variation of 
heat transfer coefficient at the walls for different 
values of the physical parameters of interest. We 
observe that nature of the heat transfer is in 
oscillatory behavior, which may be due to 



peristalsis. From figures 19, 22 and 23, the 
magnitude of heat transfer coefficient increases by 
increasing Da, Br, E h E 2 and E 3 while from Figures 
17, 18, 20, 21 and 22, it decreases by increasing 
\,p, /Wand m. 

TRAPPING PHENOMENA 

The effect of slip parameter on the streamline 
pattern is shown in Fig. 24. We observe that the 
size of the trapping bolus increases with increasing 
slip parameter. The streamlines for Jeffrey 
parameter \ are shown in Fig. 25. It is observed 
that the size of trapping bolus decreases with 
increasing \. The streamlines for uniform and non- 
uniform channels are shown in Fig. 26. It is found 
that the size of the trapped bolus is large in the left 
hand of the convergent channel while it has 
opposite behavior for divergent channel. Further, 
the size of bolus is symmetric for uniform channel. 
From Fig. 27 it can be seen that the volume of the 
trapped bolus decreases with increase of M. From 
Fig. 28 it is clear that the trapped bolus increases 
in size as Da increases. 
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Fig.2. The variation of u with y for different values of X y for fixed 
£,=1, £2=0.5, £3=0.5, £ =0. 1, P =0.0, M=2, Da=1, m=0. 



Fig.4. The variation of u with y for different values of /? for fixed 
£,=1, £2=1, £3=0.5, £=0.1, M=2, Da=2, m=0, /I, =1. 




Fig. 3. The variation of u with y for different values of Da for fixed 
£,=0.5, £2=0.5, £3=0.1, £=0.1,/? =0.2, M=2, m=0.1, ^=1. 
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Fig. 5. The variation of u with y for different values of £ for fixed 
£,=2, £2=0.7, £3=0.1, M=3, f3 =0.2, Da=2, m=0.1, A, =1- 
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Fig. 6. The variation of u with y for different values of m for fixed 
£,=0.8, £2=0.5, £3=0.5, £=0.1, jB =0.2, M=3, Da=2, \ =0.5. 
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Fig. 7. The variation of u with y for different values of /Wfor fixed 
£,=2, £2=0.7, £,=0.1, £=0.2, f3 =0.2, Da=2, m=0.1, /I, =1. 




Fig. 8. The variation of u with y for different values of £ ( , £ 2 and £ 3 
for fixed £=0.1 , f3 =0.2, M=2, Da=2, m=0, X y -\. 




Results of Sriniras et al.[20] 
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Fig. 9. The variation of 6 with y for different values of A t for fixed 
Br=3, £,=0.8, £2=0.5, £3=0.2, £=0A,f3 =0.1, M=2, Da=\, m=0.2. 
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Fig. 1 0. The variation of 6 with y for different values of Da for fixed 
Br=4, E,=1, £2=0.5, £ 3 =0.2, £=0.1, M=2, m=0.2, /? =0.1, X x =1. 



Fig.1 2. The variation of 6 with y for different values of /? for fixed 
Sr=3, £,=0.7, £2=0.5, £3=0.2, £=0.2, M=3, Da=2, m=0.1, A =0.5. 
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Fig. 11 The variation of 6 with y for different values of Mfor fixed 
Br=2, £,=1, £2=0.8, £3=0.2, £=0.1, Da=0.5, m=0.2, /? =0.1, \ =0.5. 



Fig. 13. The variation of 6 with y for different values of £ for fixed 
Br=3, £,=0.8, £2=0.5, £,=0.3, A t =1, /? =0.1, M=2, Da=1, m=0.2 
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Fig. 14. The variation of 6 with y for different values of m for fixed 
B/=4, £,=0.8, £2=0.5, £3=0.2, £=0.1, M=2, Da=1, /? =0.1, X x =0.5. 




Fig.15.The variation of 6 with y for different values of Br for fixed 
£,=0.8, £2=0.6,£ 3 =0.2, =1 , /? =0.1 , M=3, Da=\ , m=0.1 , £ =0.1 
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Fig. 16. The variation of 6 with y for different values of £,, E 2 and 
£3 for fixed Br=3, £=0.1, f3 =0.1, M=2, Da=2, m=0.1, X y = 1 . 
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Fig. 1 7. The Coefficient of heat transfer Z with x for different values 
of \ for fixed £,=0.5, £2=0.4, £3=0.1, £=0.1, M=3, B/=2, Da=0.2, 
m=0.1, P =0.1. 
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Fig. 18. The variation of Z with x for different values of /? for fixed 

£,=0.5, £2=0.4, £3=0.1, £=0.1, M=3, Br=2, Da=0.2, m=0.1, 2, 
=0.4 




Fig. 19. The variation of Z with x for different values of Da for fixed 
£,=1.2, £?=0.1, £3=0.1, £=0.1, M=5, Br=2, m=0.1, /? =0.1, ^=0.2. 




Fig. 20. The variation of Z with x for different values of M for fixed 
£,=0.5, £2=0.4, £3=0.1, £=0.1, B/=2, Da=0.2, m=0.1, /? =0.1, A, 
=0.2. 




Fig. 21 . The variation of Z with x for different values of m for fixed 
£,=0.5,£2=0.4,£3=0.1 , £ =0.1 ,B/=3,Da=0.2,/W=0.1 , /? =0.1 , Aj =0.4. 




Fig. 22. The variation of Z with x for different values of Br for fixed 
£,=0.8,£2=0.4,£3=0.2, £=0.1 ,M=3,Da=0.2,m=0.2, /? =0.1 , ^ =0.4. 
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Fig. 23. The variation of Z with x for different values of £,, £?, £3 
for fixed M=3, Br=3, Da=0.2, m=0.2, /? =0.1 , =0.4. 
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Fig. 28. Streamlines for (a) Da=0.01, (b) Da=0.1, (c) Da=°° with 
£,=0.6, E 2 =0.4, E 3 =0.1 , £ =0.1 , jB =0.1 , m=0.1 , M=4, \ =1 , t=0.1 . 



CONCLUSION 

In this paper, we investigated the Influence of 
slip conditions, wall properties and heat transfer 
on MHD Peristaltic transport of a Jeffrey fluid in a 
non-uniform porous channel under the 
assumptions of long wavelength and low 
Reynolds number. The analytical expressions 
are obtained for the velocity, stream function and 
temperature. The main observations of this study 
are as follows: 

• The velocity profile decreases with an 
increase in Jeffrey parameter \. 

• The temperature decreases with an 
increase in 4,/?, M and E 3 while it 
increases with increase in Da, e, m, E-, 
and E 2 . 

• As expected, the coefficient of heat 
transfer is oscillatory in nature. 

• The size of trapped bolus is smaller in 
Jeffrey fluid when compared with that of 

Newtonian fluid (\ = 0) . 

• As the Jeffrey parameter/?, -^0, the 
results deduced are found to be in 
agreement with the corresponding ones of 
Srinivas et al. [20]. 
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